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The liquid-gas phase transition is analyzed from the topologic properties of the event distribution 
in the obervables space. A multi-canonical formalism allows to directly relate the standard phase 
transition with neutral particles to the case where the non saturating Coulomb interaction is present, 
and to interpret the effect of Coulomb as a deformation of the probability distributions and a rotation 
of the order parameter. This formalism is applied to a statistical multifragmentation model and 
consequences for the nuclear multifragmentation phase transitions are drawn. 

PACS numbers: 24.60.-k,21.10.Sf,68.18. Jk,03.65.Vf 



Since the first pioneering multifragmentation experi- 
ments 1], the break-up of a hot nucleus above about 
3A.MeV excitation energy has been tentatively associ- 
ated to a transition from liquid to gas. A first order 
phase transition can be unambiguously defined in a fi- 
nite system as an abnormal curvature of a thermody- 
namic potential e.g. a convexity of the entropy 0- When 
this anomaly affects the energy axis, the heat capacity 
becomes negative. The first experimental signatures of 
negative heat capacities have been recently reported for 
atomic clusters 0,0] and nuclei [HE|- The results on Na- 
clusters have been related to the melting of Aa-matter^. 
For nuclei the link with bulk properties is problematic 
since the long range Coulomb interaction forbids the def- 
inition of a thermodynamic limit. Nuclear matter is a 
theoretical construction assuming the proton charge to 
be zero and one may worry about the connection between 
the properties of charged nuclei and the thermodynamics 
of uncharged nuclear matter. In this letter, we propose 
a systematic way to perform such a connection. 

Let us develop the link between charged and un- 
charged systems. The long range Coulomb interaction 
can be written as E c = q 2 V c where q represents the 
proton charge (in electron charge unit q = 1) while 
v c = J2 P<P ' a l \ r P - V I witn a = e 2 /4-7re . The total 
energy of the system E can thus be splitted into two 
terms E to t = E n + E c = E n + q 2 V c . If the system is 
in contact with a heat bath (canonical ensemble) we can 
sort the events as a function of both E n and V c : 

Pp(E n , V e ) = -^W(E n , V c )e-^-^ v '. (1) 

The canonical total energy distribution is an integra- 
tion of the generalized distribution QJ: pp(E tot ) — 
J dV c pp(E tot — q V c ,V c ). The microcanonical entropy, 
Stot(Etot) = \ogW to t(Etot), is also directly related to 
Pp(E tot ): S to t{Etot) = log P/3 (-Etot) +PE tot +logZ /3 . Eq. 
{JU naturally leads to the definition of a generalized en- 
tropy S(E n ,V c ) — logW(E n ,Vc)- The standard entropy 
can be obtained by a projection of the density of states 



on the total energy E to t — E n + q 2 V c axis 

W tot {E tot ) = J dV c W{E tot - q 2 V c , V c ) 

One can also notice that the projection of W(E n , V c ) on 
the E n axis provides the entropy of the uncharged system 

S n {E n ) = log W n {E n ) = log j dV c W{E ni V c ) 

Since we have introduced two observables {E n , V c ), it is 
natural to extend the canonical ensemble by introducing 
two Lagrange multipliers f3 n , f3 c leading to the distribu- 
tion of events in the two-dimensional observables space 

Pf 3 nf3c (E n , V c ) = -^—W{E n , Vc )e-^ E "-^ 2y c (2) 

where Zp n p a is the partition sum of this multi-canonical 
ensemble corresponding to the independent observation 
of the two components E n , V c of the energy. In the cases 
/3 n — Pc eq. J21 reduces to the usual canonical distribution 
of interacting charged particles while (3 C — corre- 
sponds to systems not affected by the Coulomb force. All 
intermediate values < (3 C < (3 n give rise to interpolat- 
ing ensembles between the charged and the uncharged 
case, in the same way as the gaussian ensemble gives 
a continuous interpolation between the microcanonical 
and the canonical ensemble. The multicanonical ensem- 
ble can also be considered as the canonical ensemble at 
temperature (3~ l of particles charged with an effective 
charge q 2 /q 2 = f3 c /(3 n - In fact we can introduce an effec- 
tive energy as E e ft — E n + q 2 V c . Then the distribution 
of the effective energy is a canonical distribution with 
a Boltzmann factor e~@ nEeff . Considering (3 C = (or 
(3 C — > faster then the increase of the charge when the 
volume of the system is increased), the multicanonical 
ensemble is a way to study the thermodynamic limit. 

Let us now focus on the question of (first order) phase 
transitions by studying the convex region of the en- 
tropy S{E C) E n ), Eq.© shows that the curvature ma- 
trix of the entropy and of the probability distribution 
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logpp n p c (E c , V n ) are the same for every set of Lagrange 
multipliers P n ,(3 c - This implies that the whole thermo- 
dynamics of the uncharged system can be (at least in 
principle) completely reconstructed from the only knowl- 
edge of the distribution probability of the charged system 
at an arbitrary temperature and pressure. For example, 
the uncharged system entropy can be deduced within a 
constant from the statistical distribution of the charged 
system by projecting out the Coulomb energy after a 
proper Boltzmann weighting 

S n (E n ) = log ( [ dV B ^ v 'p Me {E n , V c )) - (3 n E n . 



(3) 

The only limitation can be a practical one for experi- 
ments or simulations: to accumulate enough events at 
every location. 

In order to illustrate these general ideas and make 
some connections to the nuclear multifragmentation ex- 
periments, we have made some calculations in the multi- 
canonical ensemble with a statistical multifragmentation 
model ||. The model describes a multifragmentation 
event as an ensemble of spherical, non overlapping frag- 
ments interacting through their mutual Coulomb repul- 
sion. The masses of the fragments are parametrized from 
a charged liquid drop model while their excitation en- 
ergies are taken from a Fermi gas distribution with a 
high energy cut-off. The fragment translational degrees 
of freedom are treated classically. Each event is char- 
acterized by its total mass and charge number (i.e. the 
size of the disassembling nuclear source) , its total energy 
that can be decomposed in a Coulomb E c part that in- 
cludes also Coulomb effects on the binding energy, and 
a non Coulomb one that corresponds to the nuclear in- 
teraction inside fragments plus the translational energy. 
Each event can also be associated to a spatial extension, 
given by the box ('freeze-out' volume) that contains all 
the fragments. In an open system as a multifragment- 
ing nucleus, this volume is not fixed but can fluctuate 
from one event to the other. If the average volume at 
freeze-out is finite, we can treat the volume as an extra 
observable known in average taking this constraint into 
account as an additional Lagrange multiplier A ~ f3P 
where (3 = (3 n is the inverse t emp erature and P has the 
dimension of a pressure 0, llCj . The multi-canonical 
probability J2J is then evaluated with a Metropolis tech- 
nique @. 

The resulting distribution of events can then be plotted 
in the (E n ,V c ) plane or equivalently in the (E — E n + 
V c , V c ) plane. When the Coulomb interaction is at play 
(q = 1) the total energy is simply the E axis. When 
the Coulomb energy is not effective (q = 0) the relevant 
energy axis is the E — V c one i.e. the second bisector. 

In the isobar ensemble we are considering, energy is ex- 
pected to be an order parameter for the ordinary liquid- 
gas phase transition (in the uncharged system), imply- 



ing that, for all pressures below the critical one and for 
all temperatures in the transition region, we expect the 
event cloud to separate into two components (phases) 
along the energy axis Q ■ 



FIG. 1: Event distribution in the Coulomb energy ver- 
sus total energy plane and relative projections with A = 
3 • lCP 4 /m -3 for a system of total mass number A = 50 
and atomic number Z — 23 in the isobar multicanonical 
multifragmentation model. Levels of grey and full lines: 
fin 1 = 3.7MeV, (3 C = 0; contour lines and dashed lines: 
Pn 1 = Pc 1 = S.lMeV. 

Figure 1 shows the distribution of events for a system of 
50 particles at a subcritical pressure j3P = 3 • 10 _4 /m -3 
in the Coulomb energy V c and total energy E (calculated 
respect to the ground state of the source) plane. The 
charged case (3 C = (5 n and the uncharged one /3 C = are 
displayed together with their projections on the two axes. 

Let us look at the uncharged case first. The topology 
of events is characteristic of a first order phase transi- 
tion, with an accumulation of events at low excitation 
energy and high Coulomb energy (compact configura- 
tions or " liquid" phase) and an accumulation at high en- 
ergy and low Coulomb (rarefied configurations or "gas" 
phase). The distance between the two peaks along the 
nuclear energy axis E — V c , measures the latent heat of 
the transition while the region of inverted curvature is 
related to a negative heat capacity in the microcanon- 
ical ensemble. The representation in the (E, V c ) plane 
reveals that the liquid peak is constituted by a series 
of structures. They can be attributed to different chan- 
nel openings associated to similar excitation energies but 
corresponding to different proton partitions leading to 
different V c . They constitute the microscopic origin of 
the global collective phase change that in the bulk limit 
will become the ordinary liquid-gas phase transition. 



FIG. 2: Event distribution in the Coulomb energy versus 
total energy plane with A = 5 ■ 10 -3 /m -3 for a system of 
total mass number A — 50 and atomic number Z — 23 in 
the isobar multicanonical multifragmentation model. Levels 
of grey: fi^ 1 = 7.5MeV, (3 C = 0; contour lines in the left 
panel: fi^ 1 = fi^ 1 = 6.7MeF; contour lines in the right panel: 
distribution with Coulomb reweighted to get the uncharged 



The event distribution of the charged system in Figure 
1 shows globally the same structure as in the uncharged 
case, with two important differences. First, to have the 
same height in the two relative maxima the temperature 
must be lower. The decrease of the transition temper- 
ature is due to the repulsive character of the Coulomb 
interaction [ToT ]. Second, even if the phase transition is 
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still clearly visible in the bidimcnsional representation, 
this is not true any more for the projection over the total 
energy axis (E). Indeed, the fact that the distance be- 
tween the peaks (latent heat) is comparable to the width 
of each peak implies that the bimodality is hardly per- 
ceptible in the energy direction. The usual interpreta- 
tion is then the lowering of the critical point due to the 
Coulomb interaction. The complete bidimcnsional infor- 
mation of Figure 1 demonstrates that this conclusion is 
not completely correct. Indeed while the event distribu- 
tion is almost unchanged, the projection axis is different 
in the charged case and the best order parameter (i.e. 
the direction that separates at best the two phases) is 
now almost perpendicular to the energy direction. Intro- 
ducing a better observable playing the role of the order 
parameter, O = aE + bV c , would thus restore the over- 
all picture of the first order phase transition even in the 
charged system. 

The introduction of the charge is just a rewcighting 
in the occupation probabilities of the very same en- 
tropy surface. Indeed, if we consider a charged system 
at the temperature described by the distribution 
Pft = PPn=0i,0c=0n the distribution p^ = pp n=02 ^ c = o 
at the temperature /J^ 1 in the case where the Coulomb 
interaction is not active is given (within a normalization) 
by 

Pi (E n , V c ) cx p% x {En, V c ) e (ft-fc)s»+fc^ e (4) 

An example of the reconstruction of an uncharged 
distribution, starting from a calculation including the 
Coulomb interaction is shown in figure 2. The perfect 
reproduction of the distribution demonstrates that, even 
if the inclusion of Coulomb leads the system to occupy 
the phase space in a different way, the density of states 
W itself does not change and the same information about 
its convexity anomalies and phase transitions properties 
can be inferred from the two calculations. 



FIG. 3: Event distribution in the Coulomb energy versus 
total energy plane with A = 1.26 • 10 /m for a system of 
total mass number A = 200 and atomic number Z = 82 in 
the isobar multicanonical multifragmentation model. Levels 
of grey: f}^ 1 = 5.6MeV, /3 C = 0; contour lines (from top to 
bottom): /3' 1 = = 3., 3.7, 4.4, 4.8MeV. Black contour 
lines: = 5.MeV and partially screened charge q — 0.14. 

From a practical point of view, this reconstruction is 
possible only for the energy states where the population 
sample is statistically significant; if the probability dis- 
tributions are concentrated in very different locations in 
the observables space the size of the sample needed will 
become increasingly (and soon prohibitively) high. This 
more complicated situation happens when the average 
partitions are very different in the charged and uncharged 



case. This is indeed the case if we consider heavily 
charged systems, for which the Coulomb distortions get 
more appreciable. Figure 3 shows the event distribution 
for a heavy system composed of A = 200 particles. When 
the Coulomb interaction is neglected the situation gets 
closer to the expected behavior in the bulk: the channel 
openings are not recognizeable any more and only the col- 
lective liquid-like to gas- like state change survives, with 
a minimum between the two phases getting deeper. If 
Coulomb is switched on, the event topology drastically 
changes. A much smaller portion of phase space is popu- 
lated by any calculation at a single temperature, and this 
stays true if the pressure is changed. The event distribu- 
tion turns around the first order phase transition without 
diving into the coexistence region. The external edge of 
the anomalous curvature region is only touched by the 
calculations at the lowest temperature, pointing to an 
order parameter almost perfectly perpendicular to the 
energy direction. This residual bimodality corresponds 
to the opening of asymmetric fission, and is the only rest 
of the liquid-gas phase transition in heavily charged sys- 
tems. Since the average transformations in the (E, V c ) 
plane are so far away, the reconstruction of the neutral 
matter phase diagram requires a prohibitively high statis- 
tics and cannot be performed for too heavy systems. 

In conclusion, in this letter we have introduced a multi- 
canonical formalism that allows a direct mapping be- 
tween a system composed of particles with or without an 
electric charge. This general framework allows to study 
the effect of a long range non saturating force as the 
Coulomb interaction on a first order phase transition. 
Such an interaction is seen to deform the distribution of 
events in the space of observables but its strongest effect 
is to rotate the energy axis away from the order parame- 
ter. This effect is responsible for the apparent reduction 
or suppression of the coexistence region. If the system is 
not too heavily charged, the nuclear energy remains well 
correlated with the order parameter and the transition 
phenomenology is not strongly affected by the presence 
of the Coulomb field. The Coulomb simply manifests 
itself as a lowering of the critical parameters which is 
normal since the binding energy is reduced. The convex- 
ity anomalies of the entropy can be traced back from the 
event distribution independently of the Coulomb interac- 
tion and the equations of state of the uncharged system 
can be obtained from a charged population sample if the 
Coulomb energy is measured on an event-by-event ba- 
sis and the relevant phase space is sampled with enough 
statistics. 

For heavily charged systems the considered statisti- 
cal model predicts that, at low temperature, the events 
strongly differ from the uncharged partitions. Therefore 
the two thermodynamics can only be related at high en- 
ergy when the gas phase is reached. Moreover in the 
charged system the transition from the liquid to gas 
goes around the coexistence region and so corresponds 
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to a cross over and not to a first order phase transition. 
However, the size at which this phenomenon appears is 
model dependent and more investigations are needed to 
get quantitative predictions. 

Finally we would like to stress that all the interpolating 
ensembles with < (3 C < j3 n (see black contour lines in 
figure 13 may also be interesting to model experimental 
situations in which the different parts of the energy are 
not equilibrated. This may happen in actual reactions 
since the relaxation time depends critically on the range 
of the force For example, many scenarii of mul- 

tifragmentation involve spinodal instabilities which are 
weakly sensitive to the Coulomb force Being very 

chaotic, spinodal decomposition may then well populate 
a complete phase space for which E n is the important 
physical quantity and not E to t ■ This would correspond 
to a (3 C w case. Such a situation might be suggested by 
the insensitivity to the total system charge of the frag- 
mentation pattern observed in multifragmcntation data 
from central reactions and predicted in spinodal decom- 
position simulations However, this question requires 
a strong theoretical and experimental effort which goes 
beyond the scope of the present article. 
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